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Abstract 



Relativistic heavy ions are copious sources of virtual photons. The large 
photon flux gives rise to a substantial photonuclear interaction probability 
at impact parameters where no hadronic interactions can occur. Multiple 
photonuclear interactions in a single collision are possible. In this letter, we 
use mutual Coulomb excitation of both nuclei g for moderate impact 

parameter collisions. We calculate the cross section for coherent vector meson 
production accompanied by mutual excitation, and show that the median 
impact parameter is much smaller than for untagged production. The vector 
meson rapidity and transverse momentum distribution are very different from 
untagged exclusive vector meson production. 
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The large charge (Z) of heavy nuclei gives rise to strong electromagnetic fields. For rela- 
tivistic nuclei, these fields may be treated as an almost-real virtual photon beam, following 
the Weizsacker- Williams method. In relativistic heavy ion collisions, the photon field of one 
nucleus can produce a photonuclear interaction in the other. Many types of photonuclear 
interactions are possible: nuclear excitation, coherent or incoherent vector meson produc- 
tion, and other incoherent photonuclear interactions When the impact parameter b is 
significantly larger than twice the nuclear radius Ra, hadronic interactions are not possible, 
and photonuclear reactions can be cleanly detected. When b > 2Ra, the absolute proba- 
bility for many of the photonuclear reactions can be substantial, and consequently, multiple 
reactions are possible. For example, each ion can emit a photon, exciting the other nucleus, 
leading to mutual excitation [0] |^. The excited ions typically decay by emission of one 
or more neutrons which move in the longitudinal direction with approximately the same 
momentum as the beam. This process has a distinctive event signature. 

In this Letter, we consider vector meson production accompanied by mutual Coulomb 
dissociation, as shown in Figure 1. More generally, we discuss the use of mutual dissociation 
as a 'tag' of impact parameter. The vector meson calculations follow the model for exclusive 
production in Ref. 0. We will show that vector meson production accompanied by mutual 
excitation has a different impact parameter distribution than exclusive production, and this 
alters both the rapidity and transverse momentum distributions. 

The diagram in Figure |^, corresponding to exchange of 3 photons, is dominant for vector 
meson production in coincidence with Coulomb breakup. For factorization to hold, these 
photons must be emitted independently. This was demonstrated by S. N. Gupta 0. 

Calculations in Ref. [0 have shown that the probability of a photon leaving the emitting 
nucleus in an excited state is small for the case of heavy ion two-photon particle production. 
It is even less likely that the scattering represeneted by the dashed line in Fig. |I| will leave 
the emitting nucleus in an excited state. One important difference between two-photon 
interactions and coherent vector meson production is that the meson scatters from neutrons 
as well as protons. In a GDR the protons and neutrons oscillate against each other, so a force 
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with equal coupling to protons and neutrons should not excite a GDR. The equi-coupling 
to neutrons and protons also holds for the C = + mesons (mostly the fi and 02) which can 
interact with a photon to produce vector mesons. The elastic scattering should thus not 
collectively excite the nucleus, and Fig. will be the dominant process. 

These theoretical arguments are supported by data. The STAR collaboration has ob- 
served p production with and without nuclear excitation, in gold-gold collisions at an energy 
of y/SNN = 130 GeV per nucleon. The p data with excitation was collected with a trigger 
based on mutual excitation, while the exclusive p sample was collected with a low-multiplicity 
charged particle trigger ||^. Both p spectra are similarly peaked for px < 100 MeV/c. 

The Weizsacker- Williams photon spectrum at a perpendicular distance b from the center 
of the emitting nucleus is [p|-[TT[|: 



where a ~ 1/137 is the electromagnetic coupling constant, k is the photon energy, k± the 
photon transverse momentum with respect to the direction of nuclear motion. Here, Ji is a 
Bessel function, and 7 the Lorentz boost in the target frame. For a collider, 7 is related to 
the lab- frame boost 7cm by 7 = 27^^ — 1. For heavy nuclei, the nuclear charge form factor 
F can be analytically modelled by the convolution of a hard sphere with a Yukawa potential 
of range 0.7 fm |^. 

The probability of having a mutual Coulomb excitation in a collision with impact param- 
eter b is calculated along the same lines as in 0. For low energy photons {k <~30 MeV in 
the target frame), the dominant photonuclear reaction is excitation of the target to a Giant 
Dipole Resonance (GDR). This collective excitation usually decays by single neutron emis- 
sion, with multiple neutron emission also possible ||12|. Higher energy photons can excite 
the nucleus to higher collective modes or excite individual nucleons (i.e. via 7p A). More 
energetic interactions generally lead to multiple neutron emission, sometimes accompanied 
by TT emission. 

This paper will treat two general Coulomb excitation leading to the emission 
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of any number of neutrons (Xn) and an excitation to a GDR followed by the emission of 
exactly one neutron (In). The lowest-order probability for an excitation to any state which 
emits one or more neutrons (Xn) is 

Pkxnjib)=Jdk^^^a,^^^,{k). (2) 

Here, the superscript 1 shows that this is the lowest order probability. The photoexcitation 
cross section cr^^_^^* (k) is is determined by measurements at a wide range of energies . 

At small impact parameters, Pc{b) can exceed 1 and so cannot be interpreted as a 
probability. Instead, it corresponds to the mean number of excitations. The excitation 
probability may be determined by a unitarization procedure. The probability of having 
exactly N excitations follows a Poisson distribution 

{P\b)f exp{-P\b)) 
PN[b) — . (3) 

The probability of having at least one Coulomb excitation is then Pc{Xn){b) = 1 — 
exp(-P^(x„)W). 

The probability of excitation followed by single neutron emission may be similarly de- 
termined. Here, the lowest order probability is determined as in Eq. (|^), except that the 
photon energy integration is truncated at the maximum GDR energy, and data on photoe- 
mission of single neutrons is used, avoiding uncertainties about the GDR branching ratios. 
For single neutron emission, there must be a single In excitation, unaccompanied by higher 
excitations. So, Pc(in)(&) = ^c(i„)(&) exp (-P(i(6)). 

In mutual dissociation, each individual breakup occurs independently . The probabil- 
ity is thus the square of the individual breakup probabilities, i.e. Pc(XnXn){b) = {Pc{Xn){b)y 
and Pc(inin){b) = (Pc(i„)(6))^. Mutual dissociation is an experimental signature for low-6 
events [§. When a nucleus breaks up, the neutron momentum is largely unchanged, and so 
the neutrons are detectable in zero degree calorimeters downstream of the interaction point 
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Vector mesons are produced in peripheral heavy-ion interactions at ultra-relativistic en- 
ergies through coherent interactions between the electromagnetic and nuclear fields 



A photon from the field of one nucleus fiuctuates to a quark-antiquark pair and scatters 
elastically from the other nucleus, emerging as a vector meson. The cross section is sensitive 
to the vector meson-nucleon interaction cross section. The photon energy k is related to 
the final state meson rapidity, y = l/2ln(2k/Mv), where My is the vector meson mass. 
Since the photon energy spectrum depends on b, the meson rapidity distribution da/dy also 
varies with h. A vector meson can be produced at either ion and thus the two ions act as a 
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The probability for vector meson production Pv{b) is similar to that for nuclear excita- 
tion, except that Pv{b) is small, so there is no need for a unitarization process: 

Pv{b) is the probability for producing a vector meson at one nucleus. The total probability 
is obtained by doubling Py to take into account that either nucleus can emit a photon. This 
factor of 2 was not present for dissociation because Eq. ^ and our subseqent expressions for 
Pc(Xn)(yb) and Pc{in){b) were for excitation of a specific nucleus. The natural width of the 
p° is taken into account in the integration over k 

Assuming that the sub-reactions are independent, the cross section to produce a vector 
meson accompanied by mutual dissociation is 

a{AA A*A*V) = 2jd% Pv{b)PxnXn{b) exp{-PH{b)). 

Pnib) is the mean number of projectile nucleons that interact at least once: 

Pnib) = J d^fTAif- b) (1 - expi-aNNTeim • (6) 

The nuclear thickness function, T{r) is calculated from the nuclear density distribution and 
using for the total nucleon-nucleon cross section a^^ = 52 mb (88 mb) at a center of mass 
energy of 200 GeV (5.5 TeV) per nucleon pair The factor exp (— Pj/(6)) in Eq. ^thus 

ensures that the reaction is unaccompanied by hadronic interactions. For a solid-sphere 
nucleus model, the hadronic interaction probability is 1 for 6 < 2Ra and is zero otherwise. 
Here, we calculate the interaction probability from a Glauber model. 
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Table I gives the production cross sections and median impact parameters bm for the 
different tags, as calculated from Eq. ^. The XnXn and Inln cross sections are about 1/10 
and 1/100 of the untagged cross sections, respectively. 

Fig. 1^ compares the probability of p production as a function of impact parameter for 
XnXn and Inln excitation and also without requiring nuclear excitation for (a) gold-gold 
collisions at a center of mass energy ^/snn = 200 GeV per nucleon, as are found at the 
Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory, and (b) lead- 
lead collisions at ^/SNN = 5.5 TeV per nucleon as are planned at the Large Hadron Collider 
(LHC) at CERN. These curves were obtained by evaluating the integrand of Eq. (|^) at 
different b. The b distributions are very different for tagged and untagged p production; 
this is reflected in the vastly different bm in Tables 1 and 2. The Inln and XnXn spectra 
are closer, except that XnXn is more strongly peaked for b < 20/m, likely reflecting the 
increased phase space for high-energy excitations there. This difference is reflected by a 
~ 10% difference in bm- With nuclear breakup, bm is almost independent of the final state 
vector meson. 

Figure ^ shows the rapidity distribution da/dy for p and J lip production at RHIC and 
the LHC. Spectra for XnXn and Inln breakup are shown, along with the untagged da/dy. 
The da/dy are symmetric around ?/ = because either nucleus can emit the photon. Since 
the photon spectrum falls as 1/k and y = 1/2 In (2fc/My), the vector meson is usually in 
the hemisphere that the photon came from, with large \y\ corresponding to low photon 
energies. The tagged spectra have a smaller \y\ than the untagged distribution, with a 
small difference between the XnXn and Inln calculations. Breakup preferentially selects 
collisions with smaller impact parameters, with higher median photon energies, and so the 
vector mesons are produced closer to mid-rapidity. 

The form factor for a nucleus in an excited (i.e. GDR) state may be different from 
ground state nuclei. This could conceivably affect the vector meson production, particularly 
the meson pt spectra. However, this should at most be a small effect and we neglect it here. 

The meson pt spectrum also depends on b. For production at a single source (target 
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nucleus), the meson px is the sum of the photon and scattering px-, which is largely inde- 
pendent of h. The photon px comes from the equivalent photon approximation, (2), while 
the pt from the coherent scattering depends directly on the form factor, Eq. (2) of Ref. 
||16|| . The overall meson pt spectrum is affected by interference from the two production 
sources (ions). The two amplitudes add with a 6— dependent phase factor exp {ipT ■ b). At 
mid-rapidity the two amplitudes have the same magnitude but opposite sign, because of the 
negative parity of the vector meson, and 

a{pT, h) = (Ti {pt, h) [1 - cos{pt ■ b)] (7) 

— * 

where cri{pT, b) is the cross section for emission from a single source. Of course, b is unknown, 
so Eq. d^) must be integrated over b. This integration washes out the interference term in 
Eq. (1^ except for px < h/{b), where the cross section is reduced. As (b) decreases, the 
interference extends to higher and higher pt- In this way, tagging affects the pt spectrum. 

Figure || compares the pt distributions, (PN/dp^ for p and J/ip production at RHIC and 
the LHC. The solid line shows the untagged spectra, while the dashed and dotted lines are 
for XnXn and Inln, respectively. The differences are moderate at RHIC and large at the 
LHC, where the exclusive vector mesons can have much larger bm than the tagged sample. 
The different distributions can be used to study interference under different conditions. 

Although we have focused on 6— tagging vector meson production, this technique should 
also be useful for studying two-photon interactions at heavy ion colliders. Here, e^e~ pro- 
duction is of special interest. At b smaller than the electron Compton wavelength, Ac = 386 
fm, the fields are very strong, and multiple pairs production is enhanced over single pair 
production Mutual excitation could be used to select events with b < Ac, and look for 
enhanced multiple pair production. 

In conclusion, we have calculated the total cross sections and rapidity and transverse 
momentum distributions for vector mesons accompanied by nuclear breakup. The presence 
of nuclear breakup is an effective tag for events with smaller average b. These events will 
have different rapidity and pt distributions from the un-tagged events, and can be used to 
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explore the effects of different pfioton spectra and b distributions. 

Tfiis work was supported by tfie U.S. Department of Energy under Contracts No. DE- 
AC-03076SF00098 and DE-AC02-98CH10886, and by tlie Swedisli Researcli Council (VR). 
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TABLES 



TABLE I. Cross sections and median impact parameters bm, for production of vector mesons. 



Meson 




overall 


XnXn 




Inln 






a [mb] 


bm [fm] a 


[mb] 


bm [fm] 


a [mb] 


hm [fm] 






Gold beams at 


RHIC (7e„, = 


= 108) 








590 


46 


39 


18 


3.5 


19 


UJ 


59 


46 


3.9 


18 


0.34 


19 




39 


38 


3.1 


18 


0.27 


19 


J/V 


0.29 


23 


0.044 


17 


0.0036 


18 






Lead beams at LHC {^ycm = 


2940) 








5200 


280 
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19 


12 


22 


UJ 
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19 


19 


1.1 


22 




460 


220 


20 


19 


1.1 


22 




32 


68 


2.5 


19 


0.14 


21 



10 



FIGURES 




FIG. 1. The dominant Feynman diagrams for vector meson production with nuclear excitation. 
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FIG. 2. The probability of /9° production with (a) gold beams at RHIC and (b) lead beams 
at the LHC as a function of b, with XnXn (dashed curve) and Inln (dotted curve) and without 
nuclear excitation (solid curve). The Inln curve is multiplied by 10 to fit on the plot. 
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FIG. 3. Rapidity spectrum da/dy for (a) p production at RHIC, (b) J/tp production at RHIC, 
(c) p production at LHC and (d) J/V' production at the LHC. The solid line is the total production, 
the dashed line for XnXn, multiplied by 10, and the dotted line is Inln, multiplied by 100. 
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FIG. 4. The transverse momentum spectrum (fa/dp^ at mid-rapidity, y = 0, for (a) p pro- 
duction at RHIC, (b) J/^ production at RHIC, (c) p production at LHC and (d) J/^ production 
at the LHC. The sohd hne is the total production, the dashed hnc for XnXn and the dotted 
Une Inln. All of the curves assume that there is interference and are normahzed so that without 
interference dn/dpj^ = 1 at pr = 0. 
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